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THREE-DIMENSIONAL  ELECTRIC  FIELD  PREDICTIONS  OF  AN  IRON-COPPER  GALVANIC  COUPLE 


Based  on  completed  experimental  electric  field  >cans  and  the  correspo  ling 
finite  element  field  predictions,  it  appears  that  the  finite  element  numerical 
technique  presents  a  stronq  analytical  tool  in  calculating  the  nearfieid 
electric  intensity  distributions  about  active  microcells.  -rhese  calculations  were 
analytically  achieved  with  the  new  double  membrane  finite  e1ement  configuration 
representing  nonlinear  polarization  and  with  a  local  tangent  slope  (impedance) 
definition  dependent  on  the  local  potential  dif'erence.  The  experimental  determination 
of  the  multidimensional  current  density  structure  was  realized  with  a  newly 
developed  scanning  vibrating  electrode  technique  (SYCT).  The  'mite  element  model 
developed  in  this  paper  uses  a  priori  measured  uncoupled  polarization  curves  for 
pure  iron  and  pure  copper.  The  current  densities  and  the  electric  field  intensities 
were  calculated  in  the  X,  f,  and  Z  directions  wi’hin  specific  regions  of  the 
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'"'-electrolyte  and  on  its  boundaries.  Results  appear  to  indicate  that  first-order 
anodic  mass  loss  can  be  predicted  using  (1)  numerically  predicted  current 
density  distributions  on  the  anodic  surface  and  (2)  Faraday's  law.  The 
electric  field  correlation  established  in  this  work  for  the  three-dimensional 
current  density  components  provides  the  confidence  to  proceed  in  the  evaluation 
of  time-dependent  effects  of  electric  fields  and  multi  polarized  surfaces 
associated  with  pitting  and  crevice  corrosion. 
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THREE-DIMENSIONAL  ELECTRIC  FIELD  PREDICTIONS  OF  AN 
IRON-COPPER  GALVANIC  COUPLE 


INTRODUCTION 


It  Is  generally  agreed  that  electrogalvanlc  Interactions  of  anodic  and 
cathodic  regions  Inherent  In  real  world  physical  systems  are  very  complex 
processes.  This  statement  Is  particularly  true  In  view  of  the  highly 
Interactive  coupling  that  exists  between  multlp'e  anodes  and  cathodes 
subject  to  geometric  effects.  Inhomogeneous  electrolytic  mediums, 
polarization  effects,  etc.  Such  Interactions  can  be  characterized  n  a 
global  or  macroscopic  sense  as  a  complex  dc  electric  field  problem  usually 
governed  by  sensitive  nonlinear  boundary  conditions,  (l.e.,  polarization 
curves  for  anodes  and  cathodes  subject  to  the  boundary  conditions  Invariably 
enforced  In  all  boundary  value  problems  of  partial  differential  equations). 
Traditionally,  most  of  the  emphasis  In  electrochemical  research  programs  Is 
focused  on  empirical  progress  to  further  our  understanding  of  the  various 
mechanisms.  Often,  therefore,  galvanic  series  charts  oversimplify  the 
Interactions  between  a  variety  of  Ratals  and  alloys  and  thus  are  Inadequate 
and  misleading  In  terms  of  predicting  detailed  corrosion  behavior. 

A  quantitative  attempt  to  understand  these  complex  electrolytic 
mechanisms  was  made  some  40  years  ago  with  the  Introduction  and  description 
of  mixed-potential  theory.1  This  theory  state:  that  any  electrochemical 
reaction  can  be  divided  Into  one  or  more  oxidation  reactions  and  one  or  more 
reduction  reactions.  It  further  asserts  that  there  exists  no  net  charge 
accumulation  during  an  electrochemical  reaction.  This  Is,  of  course,  an 
expression  of  conservation  of  charge.  The  significance  of  these  statements 
Is  In  the  recognition  that  electrically  coupled  metals  In  an  electrolyte 
behave  totally  differently  rom  the  same  metals  when  they  are  electrically 
Isolated.  Hence,  the  observed  synergism  between  the  more  dominant  variables 
needs  to  be  somehow  mathematically  Integrated  Into  the  total  solution 
process . 


BACKGROUND 


The  traditional  method  of  determining  adequate  cathodic  protection  for  a 
physical  systjp  Is  generally  based  on  experience  and  extremely  crude 
guidelines.  Past  analytical  attempts  to  calculate  local  current 
distributions  In  the  anodic /cathodic  neighborhood  (surfaces)  generated 
Immediate  realizations  that  exact  mathematical  solutions  seemed 
Intractable.?  It  appeared  that  any  closed  form  mathematical  attempts 
would  be  restricted  to  simplified  geometric  electrochemical  sources. 
Various  theoretical  treatments  were  published  In  the  1950’s  by  Wagner3'4 
and  by  Waber,5"7  addressing  coplanar  electrode  configurations  with  varying 
degrees  of  polarization  effects.  These  analyses  were  generally  limited  to 
one  dimensional  efforts  In  Cartesian  coordinate  systems  with  simplified 
polarization  mechanisms.  A  mathematical  analysis  of  circular  corrosion 


cells  with  unequal  polarization  parameters8  demonstrated  the  treatment  of 
slightly  more  complex  geometry  In  the  presence  of  more  realistic 
polarization  effects. 

In  recent  years,  mathematical  models  have  been  developed  that  apply 
finite  element  numerical  techniques  to  the  prediction  of  macroscopic 
electrogalvanic  fields.  These  models  were  developed  for  performance 
predictions  of  cathodlcally  protected  structures,  particularly  In  the  marine 
and  nuclear-cooling  environments .9_1*  Although  the  models  did  not 
consider  changes  In  geometry  and  electrolyte  properties,  they  were 
considered  quite  successful  In  predicting  current  distributions  at  the 
various  anodic  and  cathodic  areas  of  large-scale  structures.  Furthermore, 
Improvements  In  defining  nonlinear,  grad lent -dependent  boundary  conditions 
at  electrode  surfaces  have  provided  significant  progress  In  addressing 
complex  polarization.15*18 

More  recent  advances  In  techniques  for  measuring  localized  currents  In 
solutions  have  Introduced  methods  for  measuring  Ionic  currents  associated 
with  corrosion  microcells.  The  new  techniques  use  vibrating  probes  to 
measure  currents  In  a  solution  with  spatial  and  current  density  resolutions 
on  the  order  of  15-20  yin  and  5  nA/cm2,  respectively.17  Thus,  It  Is 
possible  to  measure  corrosion  currents  resulting  from  Individual  local  cell 
activity  un  a  scale  close  to  the  size  of  many  mlcrostructural  features  of 
materials.  This  technique  Is  particularly  useful  for  studies  of  localized 
corrosion  phenomena,  such  as  pitting  and  crevice  corrosion,  and  for  studies 
of  '•orroslon  In  composite  materials  where  local  galvanic  effects  between 
constituents  and  Interfaces  may  be  Important.18*19 

Progress  In  the  finite  element  method  for  predicting  one -dimensional 
electrogalvanic  fields  using  data  from  a  one -dimensional  scanning  vibrating 
electrode  technique  (SVET)  was  documented  by  Kasper  and  Crowe.28  The 
model  that  was  constructed  provides  for  nonlinear  polarization  effects  and 
the  possibility  that  these  polarization  effects  will  evolve  with  time.  The 
results  demonstrated  the  steady-state  capabilities  of  the  model,  as  well  as 
the  need  to  resolve  the  true  three-dimensional  nature  of  the  local  electric 
field  density.  The  uniqueness  of  this  work  Is  that  the  electrochemical 
measurements  needed  to  verify  calculated  vector  components  of  an 
electrogalvanic  field  predicted  by  the  finite  element  method  were  made. 
These  measurements  were  taken  In  one-  and  three-dimensional  spatial 
configurations. 

The  mathematical  treatment  In  this  report  addresses  the  true  nonlinear 
aspects  of  the  polarization  of  the  anode  and  the  cathode  surface  areas.  The 
measured  electrode  polarizations  become,  In  effect,  nonlinear  Impedance 
relationships  that  are  functions  of  the  local  potential;  l.e.,  the  vector 
constitutive  equation  becomes  nonlinear. 


ELECTROCHEMICAL  MATERIALS  AND  MEASUREMENTS 


A  galvanic  couple  was  constructed  of  6.35-mm-diameter  Marz  grade  iron 
and  50.8-mm-diameter  Marz  grade  copper  obtained  from  MRC,  Inc.,  Orange,  New 
Jersey.  The  iron  anode  was  press  fit  into  a  centrally  located  hole  in  the 
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copper  cathode  to  produce  a  cylindrical ly  symmetrical  planar  galvanic 
couple,  as  shown  In  figure  1.  The  couple  was  metal lographlcally  polished 
using  standard  techniques,  with  T.O-pm  diamond  paste  providing  the  final 
finish.  Individual  specimens  for  polarization  measurements  were  produced 
from  the  same  materials,  and  all  were  thoroughly  cleaned  In  ethanol  before 
testing. 

Electrochemical  tests  were  conducted  In  open  cells  In  laboratory  air. 
All  electrochemical  measurements  were  made  using  0.6  N  NaCl  electrolytes 
prepared  with  reagent  grade  chemicals  and  distilled  water.  The  pH  was 
adjusted  to  1.75  with  HC1 .  The  measured  solution  conductivity,  o,  was 
63  mS/cm  at  pH  «  1.75. 

The  standard  potentlostatlc  polarization  response  of  the  constituent 
materials  was  measured.  The  specimens  for  these  tests  were  machined  from 
the  same  copper  and  Iron  rods  used  to  construct  the  galvanic  couple. 
Polarizations  were  performed  usl-ig  a  PAR  model  173  potentlostat  In 
combination  with  a  PAR  model  175  universal  programmer  at  a  stepping  rate  of 
0.1  mV/s.  Standard  ASTN  65-78  procedures  were  used. 

Ionic  current  density  measurements  were  performed  using  the 
three-dimensional  vibrating  probe  Instrument  shown  schematically  In 
figure  2.  This  Instrument  was  designed  and  built  by  E  and  B  Systems, 
Nashville.  Tennessee.  It  has  six  major  components:  the  microcomputer,  the 
drive  heads  and  amplifiers,  the  microelectrodes ,  the  microelectrode 
preamplifier,  and  two  constant  current  sources. 

The  IBH-PC  microcomputer  provides  Instrument  control,  data  acquisition, 
signal  processing,  and  display  and  storage  functions.  Instrument  control 
and  data  acquisition  are  accomplished  by  computer  Interface  to  a  Tekmar  Lab 
Haster  AD/DA  converter.  In  addition  to  controlling  the  motion  of  the 
vibrating  probe  and  the  measurement  currents,  the  computer  and  Interactive 
software  provide  for  calibration  of  the  system,  measurement  of  the  electrode 
Impedance,  and  data  analysis  and  output. 

The  drive  heads  consist  of  three  electromagnetic  drivers,  which  are 
connected  to  the  measuring  microelectrode  by  3-mm  mylar  tubes  arranged  In 
such  a  fashion  that  the  three  orthogonal  drives  axes  form  a  cube 
(figure  3).  The  diagonal  of  the  cube  Is  positioned  vertically.  Each  driver 
consists  of  a  modified  permanent -magnet  loudspeaker,  which  vibrates  the 
measuring  microelectrode  through  the  mylar  tubes  using  a  push-pull  motion. 
The  drivers  respond  to  sinusoidal  signals  generated  In  the  microcomputer. 
These  signals  e  e  applied  to  the  drive  heads  through  the  driver  amplifiers, 
resulting  In  a  circular  motion  of  the  probe  tip  In  a  plane  angled  at  54.7 
degrees  to  the  axis  of  the  probe.  The  relationship  between  the  X-Y-Z 
specimen  coordinate  system  and  the  A-B-C  probe  drive  coordinate  system  Is 
shown  In  figure  4.  The  frequency  response  of  the  drivers  Is  shown  In 
figure  5.  As  can  be  seen,  the  drivers  provide  adequate  displacements  for 
measurement  of  the  current  density  at  of f -resonance  frequencies.  In 
contrast,  the  piezoelectric  drivers  used  In  one -dimensional  probes  lose 
their  stroke  rapidly  when  operated  of f -resonance. 

The  microelectrode  system  consists  of  two  electrodes:  the  measuring 
electrode  and  the  reference  electrode.  The  measuring  electrode  is  a 
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platinized,  gold-plated,  tungsten  wire  electrode,  which  Is  Insulated  with 
paralene.  It  Is  mounted  In  the  junction  of  the  mylar  drive  tubes,  with  the 
tip  of  the  probe  exposed,  as  shown  In  figure  6.  The  reference  electrode  is 
a  sllver/sllver  chloride  electrode  rigidly  attached  to  the  preamp! .fler. 
The  preamplifier  contains  a  two-stage  differential  preamplifier,  which 
provides  a  gain  of  104  In  the  voltage  difference  between  the  vibrating 
measuring  electrode  and  the  reference  electrode. 

The  theoretical  sensitivity  limit  of  the  electrode  Is  determined  by  the 
electrode  Impedance  and  thermal  noise  of  the  electrode.  The  exposed  area  of 
the  electrode  tip  determines  both  Its  spatial  resolution  and  Impedance. 
Thu.,,  tip  preparation  Is  Important. 

When  the  electrode  area  Is  small,  as  required  for  fine  spatial 
resolution,  the  resistive  component  of  the  electrode  Impedance  Is  high, 
causing  an  appreciable  noise  level.  This  electrode  Impedance  depends  upon 
both  the  electrode  area  and  the  frequency  of  vibration.  The  relationship 
between  these  parameters  Is  given  by  Beddes  as  approximately21 


,  0.225  x  1014 

Z  *  - —  *  (1) 

kaf(W) 

where  Z  Is  the  electrode  Impedance,  k  Is  a  constant  that  depends  on  the  type 
and  condition  of  the  electrode,  a  Is  the  electrode  area  (In  wm2),  f  Is 
the  frequency  of  vibration,  and  a  Is  a  constant  that  relates  the  rate  at 
which  the  electrode  capacitance,  c,  decreases  with  frequency  ( 1 . e . , 
c  *  kf "•) .  The  electrode  Impedance,  measured  before  each  run,  was 
typically  1  to  1.5  x  105  Q  at  1.5  kHz  and  5  to  6  x  10s  a  at  300  Hz. 

The  ability  of  the  electrode  to  efficiently  detect  the  electric  field 
depends  upon  the  noise  level  of  the  electrcle  and  the  amplification  of  the 
signal.  Electrode  noise  arises  from  the  random  movement  of  charge  carriers 
at  the  electrode/electrolyte  Interface  as  the  electrode  establishes  and 
maintains  equilibrium  with  the  electrolyte  or  as  the  electrode  responds  to 
the  passage  of  current.  This  random  movement  of  charge  carriers  creates  a 
resistive  component  to  the  electrode  Impedance,  which  becomes  the  noise 
generator.  The  noise  voltage  Is  generated  at  all  frequencies  and  Is  given 
by  the  relation 


E  -  [4kTR(fh  -  f t ) ]0  5  .  (2) 


In  this  relation,  E  Is  the  root -mean -square  noise  voltage,  k  Is  Boltzmann’s 
constant,  T  Is  the  absolute  temperature,  R  Is  the  universal  gas  constant, 
and  fj,  and  fj  are  the  upper  and  lower  cutoff  frequencies  of  the 
measuring  system  (or  the  bandwidth). 

To  minimize  the  noise.  It  Is  necessary  to  have  as  narrow  a  bandwidth  as 
possible.  An  ultra-narrow  bandwidth  Is  obtained  during  data  reduction  using 
the  ultra -narrowband  phase  coherent  filtering  algorithm.21  This  algorithm 
allows  the  gradients,  as  well  as  the  current  density  components,  to  be 
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,  extracted  from  the  data  measured  during  rotation  of  the  electrode  in  a 

I  circle.  The  measured  voltages  from  the  rotating  electrode  are  represented 

by  a  Taylor-serles  expansion  about  the  central  point  (A0,B0)  of  the 
!  circle  of  radius  r  swept  out  during  rotation  through  angle  6  of  the  probe 

tip.  Thus, 


♦(A.B)  = 


♦<  W  *  M 


“  -  v  ♦  H 
<vv 


♦<VBo>  *  r  cos  e 


a* 

3A 


(B  -  B0) 

<VBo> 


+  r  sin  $ 


M 

38 


(VBo> 


(3) 


<VBo> 


where  (A  -  A0)  *  r  cos  6,  (B  -  B0)  *  r  sin  0,  and  the  partial  deriva¬ 
tives  are  evaluated  at  (A0,B0).  The  current  density  at  the  central 
point  (A0,B0)  has  components  3^  and  JB  91ven  By  the  convolution 
Integrals 


JA  -  “  ?  Jo*  V(6)  cos  6  d  6 

and  (4) 

JB  *  "  r  Jo*  V(6)  cos  e  d  e  • 

where  V(e)  Is  the  measured  voltage  at  rotation  angle  e  and  a  is  the 
conductivity  of  the  electrolyte  at  (A0,80).  Use  of  this  procedure 

allows  current  densities  of  5  nA/cm2  to  be  detected. 

In  operation,  the  computer  feeds  sinusoidal  driving  signals  to  the  drive 
heads.  In  sequence,  two  at  a  time,  so  that  the  probe  tip  Is  forced  to  rotate 
In  circles  In  the  three  orthogonal  planes  angled  at  54.7  degrees  to  the 
probe  axis.  The  diameter  of  the  circle  is  adjusted  as  an  experimental 

variable.  The  differential  voltage  Is  then  sampled  at  256  points  around  the 
circumference  of  each  circle.  These  data  are  averaged  over  1000  cycles  or 
more  and  then  stored  In  the  computer  for  analysis. 

During  analysis  the  microcomputer  evaluates  the  convolution  Integrals 
for  1000  or  more  cycles  of  the  probe  rotation,  adjusts  for  orientation 

relative  to  the  sample,  averages  the  current  densities,  and  calculates  the 
gradients . 

As  each  of  the  combinations  of  drives  Is  energized  in  sequence,  a 

three-dimensional  volume  of  electrolyte  Is  sampled  and  a  redundant  set  of 
the  components  of  the  current  density  vector  and  gradients  Is  obtained. 
Transf ormat Ion  of  coordinate  systems  from  the  A-B-C  to  the  X-Y-7  system  is 
performed  In  the  computer,  displayed  on  the  computer  monitor,  and  stored  in 
ASCII  files  on  the  disks  for  future  analysis. 
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Calibration  of  the  probe  is  accomplished  by  measuring  a  known  uniform 
current  density  field  in  an  electrolyte  contained  in  a  separate  calibration 
dish,  which  is  mounted  on  a  turntable.  The  field  is  rotated  through  12 
standard  positions,  30  degrees  apart.  At  each  position,  the  output  voltage 
is  measured  when  a  single  drive  is  energized.  This  technique  allows  the 
average  line  of  motion  of  the  electrode  in  the  horizontal  plane  to  be 
determined  for  that  drive,  as  well  as  the  mechanical  lag  between  the  drive 
voltage  and  the  mechanical  response  of  the  probe.  Appropriate  mathematical 
calculations  of  these  data  provide  an  exact  description  of  the  probe  motion 
and  the  measured  values  of  the  known  field  from  which  calibration  factors 
are  extracted  and  stored  for  data  reduction. 

Mechanical  positioning  of  the  probe  in  the  X-Y  plane  parallel  to  the 
specimen  surface  is  achieved  using  two  Klinger  model  MT-160  stepping 
motor-driven  translation  stages  mounted  at  right  angles.  These  translation 
stages  allow  indexing  of  the  specimen  relative  to  the  probe  in  l-pm 
increments  to  an  average  accuracy  of  2  ym.  Control  of  the  stages  is 
achieved  either  manually  or  by  computer  using  a  Klinger  model  CC-1  universal 
programmable  indexer. 

For  the  experiments  reported  here,  the  vibration  probe  tip  was 
positioned  above  the  surface  of  the  iron-copper  galvanic  couple  at  either 
650-  or  1150-ym  elevation  and  current  density  measurements  were  made.  The 
arrangement,  which  is  shown  schematically  in  figure  7,  was  designed  so  that 
the  three-dimensional  probe  measurements  were  directly  comparable  with  the 
one-dimensional  measurements  reported  earlier.^  At  each  elevition,  a 
single  line  scan  was  performed  so  that  the  probe  tip  traveled  parallel  to  a 
line  passing  through  the  center  of  the  galvanic  couple. 


EXPERIMENTAL  RESULTS 


The  anodic  and  cathodic  potentiostatic  polarization  curves  for  the 
uncoupled  ‘ ron  and  copper  (Fe-Cu)  in  0.6  M  NaCl  at  a  pH  level  of  1.75  and  a 
temperature  of  25°  C  are  shown  in  figure  8.  These  polarization 
characteristics  were  reflected  in  the  experimental  setup,  as  well  as  in  the 
impedance-control  led  boundary  conditions  of  the  finite  element  model. 


The  components  of  the  measured  ionic  current  density  vectors, 
coordinates,  as  a  function  of  position,  are  listed  in  table  1. 


polar 


Numerical  calculations  were  performed  using  the  finite  element  model 
developed  in  reference  20  and  summarized  in  figure  9.  This  finite  element 
model  contains  (1)  a  geometric  description  of  the  galvanic  couple  and  the 
electrolyte  (figure  9(a)),  (2)  the  fundamental  electrical  laws  that  govern 
the  behavior  of  the  galvanic  couple  in  partial  differential  equation  form 
(figure  9(b)),  and  (3)  the  associated  Neumann  boundary  conditions 
(figure  9(c)).  The  numerically  predicted  values  for  Jr  and  J i  are 
listed  in  tables  2  and  3  for  assumed  electrolyte  conductivities  of  42  mS/cm 
and  in  tables  4  and  5  for  conductivities  of  62.5  mS/cm.  Figure  10  shows  the 
calculated  and  measured  J^  ionic  current  density  component  as  a  function 
of  radial  position  at  650  ym  above  the  electrode  surface.  Figure  11  shows 
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Table  1.  Measured  Ionic  Current  Density 
(pH  =  1.75,  o-42  mS/cm) 


Radius 

(mm) 

Height 

(pm) 

JZ, 

(A/m2) 

Jr 

(A/m2) 

J  Total 
(A/m2) 

0 

650 

14.300 

5.430 

15.296 

1 

650 

8.420 

12.250 

14.865 

2 

650 

4.400 

11.210 

12.043 

3 

650 

2.350 

9.320 

9.612 

5 

650 

0.225 

5.350 

5.355 

7 

650 

-0.100 

2.790 

2.792 

9 

650 

-0.044 

2.030 

2.030 

n 

650 

-0.139 

1 .670 

1  .676 

13 

650 

0.007 

1 .500 

1.500 

15 

650 

-0.068 

1  .370 

1  .372 

17 

650 

-0.162 

1 .310 

1  .320 

19 

650 

0.005 

1 .180 

1 .180 

21 

650 

-0.009 

1 .090 

1  .090 

23 

650 

-0.142 

0.990 

1  .000 

25 

650 

-C.007 

0.710 

0.710 

1150 

2.960 

7.517 

1 

1150 

3.520 

4.320 

5.573 

2 

1150 

2.310 

4.480 

5.040 

3 

1150 

1  .660 

4.250 

4.563 

4 

1150 

0.858 

3.980 

4.071 

6 

1150 

0.407 

3.070 

3.097 

8 

1150 

0.111 

1 .650 

1  .654 

10 

1150 

0.608 

1  .050 

1  .213 

12 

1150 

0.061 

1.120 

1.122 

14 

1150 

0.049 

0.940 

0.941 

16 

1150 

-0.008 

0.940 

0.940 

18 

1150 

-0.050 

0.950 

0.951 

20 

1150 

0.071 

0.840 

0.843 

22 

1150 

-0.004 

0.800 

0.800 

24 

1150 

0.051 

0.770 

0.772 

0 

650 

14.300 

5.430 

15.296 

0 

1150 

6.910 

2.960 

7.517 

0 

1650 

5.290 

7.080 

8.830 

0 

2150 

4.360 

1  .850 

4.730 
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Tabic  2.  Nuaerlcal  Modeling  Results 
(Z  •  650  *m,  9  -  42  aS/ca) 


Radius 

(■a) 

Jr 

(A/*2) 

Jz_ 

(A/a2) 

3  Total 
(A/a2) 

0.7 

♦•2.414 

10.57 

10.842 

1  .786 

3.7582 

7.4759 

8.367 

- 

2.778 

3.6514 

4.8028 

6.033 

4.375 

3.3654 

2.5709 

4.235 

5.897 

2.6823 

0.93195 

2.839 

• 

7.713 

1.5649 

0.0633 

1 .566 

9.73 

0.98767 

0.1707 

1 .002 

11.29 

0.689 

0.1820 

0.7126 

13.875 

0.44 

0.1975 

0.4822 

17.125 

0.2558 

0.2232 

0.33948 

20.375 

0.1185 

0.2239 

0.2533 

23.65 

0.03638 

0.15278 

0.15705 

Table  3.  Nuaerlcal 

Modeling  Results 

(Z  -  1125  ym. 

o-42  aS/ca) 

Radius 

Jr 

Jz 

J  Total 

(«•) 

(A/a2) 

(A/a2) 

(A/a2) 

0.7 

1  .497 

7.596 

7.742 

1  . 786 

2.967 

5.989 

6.684 

2.778 

2.9182 

4.1926 

4.726 

• 

4.375 

2.894 

2.293 

3.692 

5.897 

2 . 3986 

0.9784 

2.591 

7.713 

1  .486 

0.10215 

1.4895 

9.73 

0.96688 

0.10625 

1  .002 

6 

11.29 

0  6828 

1.7137-1 

0.7039 

|  13.875 

0.436 

0.19835 

0.479 

»  17.125 

0.25559 

0.20559 

0.328 

!  20.375 

0.1195 

0.202514 

0.235 

|  23.65 

0.1641 

0  13939 

0.215 

8 
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Table  4.  Numerical  Modeling  Results 
(Z  -  650  *a,  «  -  42.5  mS/cm) 


Radius  Jr  Jz  J 

(■a)  (A/**)  (A/»2) 


0.7 

5.039 

15.99 

16.765 

1.786 

6.833 

8.2107 

10.682 

2.778 

4.8525 

2.46 

5.44 

4.375 

2.593 

0.6134 

2.664 

5.897 

1.33 

0.0568 

1 .33 

7.713 

0.8353 

0.07133 

0.83836 

9.73 

0.55485 

0.010757 

0.56518 

11.29 

0.3961 

0.10948 

0.406 

13.875 

0.2557 

0.11634 

0.2809 

17.125 

0.1507 

0.1316 

0.19959 

20.375 

0.070059 

0.1308 

0.14835 

23.65 

0.020349 

0.0885 

0.0908 

I 

! 


Table  5. 
(Z  - 


Nunerlcal  Modeling  Results 
1125  um.  a  -  62.5  «S/c«) 


Radius 

Jr 

(A/«2) 

Jz 

(A/*?) 

J  Total 
(A/m*) 

0.7 

3.0225 

10.1026 

10.545 

1.786 

5.1936 

6.406 

6.246 

2.778 

3.8138 

2.722 

4.685 

4.375 

2 . 3307 

0.813686 

2.468 

5.897 

1.323 

0.1998 

1.338 

7.713 

0.8142 

0.1831 

0.8345 

9.73 

0  5468 

0.081539 

0.5528 

11.29 

0.3931 

0.1088 

0.4078 

13.875 

0.2544 

0.12016 

0.2813 

17.125 

0.15058 

0.1226 

0.194 

20.375 

0.07069 

0.12016 

0.139 

23.65 

0.009918 

0.08258 

0.083 
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a  comparison  of  tha  measurements  at  Z  -  1150  pm  with  the  calculated 

results  at  Z  -  1125  % m.  The  agreement  between  the  Measurements  using  the 
one -dimens tonal  (1-0)  probe  (small  Insert  In  figure  11)  Is  excellent,  and 
both  sets  of  measurements  are  In  excellent  agreement  with  the  finite  element 
model  calculations.  As  can  be  seen  In  these  figures,  the  anodic  current 

density  normal  to  the  couple  Is  highest  near  the  center  of  the  Iron.  The 
normal  current  density  falls  rapidly  to  a  low  value  just  past  the 
Iron-copper  Interface,  and  there  Is  a  dispersion  In  J  as  the  elevation  above 
the  surface  Increases. 

Figure  12  shows  the  measurements  of  Jp  obtained  using  the 
three-dimensional  probe  compared  with  the  finite  element  calculations  at 

Z  -  650  iaa.  It  should  be  pointed  out  that  the  measured  components  of 
.1(JX,  Jy,  J7)  were  changed  to  a  set  of  Jp  and  J;  components  to 

compare  with  the  finite  element  results.  The  X  and  Y  components  of  J  were 
squared,  and  the  square  root  of  the  sum  of  the  squares  resulted  In  Jp. 

This  change  allowed  a  clearer  point-to-point  comparison  between  theory  and 
measurement.  (Tables  1  through  5  already  Include  this  Information.)  The 

measured  results  at  1150  tee  are  shown  compared  with  the  finite  element 

calculations  at  1125  pm  In  figure  13.  Since  the  one -dimensional  probe 
cannot  measure  these  components  of  the  current  density  vector,  there  are  no 
comparable  one -dimensional  measurements.  Again,  the  agreement  between  the 
measurements  and  the  calculations  Is  excellent. 

Figures  14  and  15  show  the  measured  and  calculated  total  current  density 
versus  radial  distance  at  Z  *  650  pm,  respectively.  The  Ionic  currents 

flowing  from  the  Iron  are  readily  discernible.  The  finite  element  values 
and  the  measured  values  from  the  three-dimensional  probe  are  plotted  in 
figures  11  through  15.  The  Iron-copper  Interface  is  located  at  a  radial 
distance  of  3.175  mm  from  the  center.  The  measured  values  from  the  second 
scan  at  1150  pm  above  the  electrode  surface,  as  well  as  the  finite  element 
results  at  1125  pm  above  the  electrode  surface,  show  excellent 
correlation.  At  both  elevations.  It  appears  that  the  current  density 
emanating  from  the  anode  surface  Is  quite  high,  with  a  fairly  evenly 
distributed  current  density  response  on  the  copper  surface.  The  total 
calculated  current  leaving  the  anodic  surface  was  +2.857612  x  10~5  A,  and 
the  total  current  arriving  over  the  copper  surface  was  -2.857670  x 
10~5  A.  This  highly  accurate  match  In  magnitude  between  the  anodic  and 
cathodic  currents  illustrates  the  accuracy  of  the  current-accounting  system 
in  the  mode  I . 


NUMERICAL  MODELING  SUMMARY 


An  extensive  discussion  of  the  numerical  procedures  for  the 
three-dimensional  finite  element  model  Illustrated  in  figure  9  has  already 
been  presented  by  Kasper  and  Crowe. 20  With  the  emphasis  here  on  the 
spatial  comparison  of  current  density  vectors  '.n  cylindrical  coordinates 
very  near  the  electrode  surfaces.  It  is  important  that  several  aspects  of 
the  mathematical  procedures  be  discussed  because  of  the  sensitivity  of  the 
current  density  vector  field  response  to  nonlinear  polarization  effects, 
choices  of  Interpolation  functions,  manner  of  discretization,  geometry, 
etc.  The  Introduction  of  the  double  nonlinear  convective  film  menbrane  to 


10 


TR  7849 


represent  the  nonlinear  polarization  of  Iron  and  copper  Is  fully  discussed 
elsewhere.15 


Yet,  some  discussion  on  the  finite  element  modeling  subtleties  Is 
essential.  The  finite  element  model  developed  for  this  Investigation  was 
constructed  from  three-dimensional  elements  using  the  COSWIC/NASTRAM* 
program.  Because  three-dimensional  finite  elements  are  more  difficult  to 
Implement,  careful  numerical  tracking  for  accuracy  and  consistency  Is 
required  to  develop  numerical  confidence.  The  results  reported  here  provide 
a  first  step  In  exercising  the  model  for  three-dimensions.  Before  their 
presentation,  however,  It  is  worthwhile  to  discuss  several  aspects  of  the 
model  and  the  resulting  computations. 

First,  the  electrolyte  was  represented  with  CHEXA2  and  CWEDGE  elements. 
The  CHEXA2  element  Is  a  three-dimensional  solid,  with  eight  vertexes  and  six 
quadrilateral  faces,  that  Is  represented  by  two  overlapping  tetrahedra.  The 
CWEDGE  element  Is  a  three-dimensional  solid  with  three  quadrl lateral  faces 
and  two  opposing  triangular  faces.  These  elements  use  linear  Interpolation 
functions  for  the  potential  between  adjacent  grid  points.  The  general 
geometric  configurations  of  the  electrolyte  and  the  anodic/cathodic  surfaces 
have  been  Illustrated  in  figure  10. 

The  model  was  bounded  by  the  following  cylindrical  coordinates:  0.0  <  R 
<  0.044  m,  0°  <  6  <  36°,  and  0.0  <  2  <  0.0095  m.  The  selection  of  the 
36-degree  wedge  angle  was  governed  mainly  by  the  choice  of  the  CWEDGE  and 
CHEXA2  elements  and  the  anodic  area  in  such  a  way  that  the  total  number  of 
three-dimensional  elements  would  be  reasonable  for  computational  purposes 
(figure  9(a)).  The  aspect  ratios  for  the  elements  appeared  reasonable, 
although  somewhat  greater  discretization  in  the  radial  direction  (especially 
in  the  anodic  area)  would  have  been  more  desirable.  On  all  surfaces  of  the 
model,  the  normal  derivative  of  the  potential  function  was  assumed  to  be 
zero  (insulated).  This  assumption  includes  the  surfaces  6=0°  and  e  - 
36°  where  the  normal  derivative  is  the  6  direction. 


The  idea  of  dealing  with  a  geometric  singularity,  such  as  a  pit  or  a 
crevice,  appeared  to  rule  out  the  possibility  of  using  any  axisymmetric 
finite  elements.  Even  though  similar  techniques,  when  applied  to 
axisymmetric  problems,  appeared  to  provide  reasonable  results,  axisymmetric 
finite  elements  are  not  adequate  in  addressing  the  geometric  changes  of  pits 
or  crevices  evolving  with  time.  Therefore,  it  was  necessary  to  introduce 
the  three-dimensional  elements  that  would  ultimately  be  required  to  address 
general  pitting  and  crevice  corrosion. 

Furthermore,  if  three-dimensional  elements  such  as  CWEDGE  and  CHEXA2 
were  used  in  model  development,  the  question  of  whether  these  elements  are 
mathematically  accurate  enough  to  capture  the  various  potential  gradient 
changes,  especially  near  the  iron  electrode,  required  attention.  The  CWEDGE 
and  CHEXA2  elements  use  linear  interpolation  functions  in  such  a  way  that 


•NASTRAN  is  an  acronym  for  the  National  Aeronautics  and  Space 
Administration  Structural  Analysis  Computer  Program.  COSMIC/NASTRAN 
indicates  that  it  is  a  Government -proprietary  NAS1RAN  code. 
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the  potential  Is  expressed  as  a  linear  algebraic  equation  In  R,  e,  and  Z. 
Another  advantage  Is  that  these  three-dimensional  elements  do  not  suffer 
from  singularity  problems  because  of  the  local  formulation  of  the 
Interpolation  function.  As  constructed,  the  model  provides  the  continuous 
variation  of  e  as  a  function  of  Z  at  R  -  n 

Finally,  It  was  desirable  to  attempt  to  control  the  current  flow  through 
the  polarization  layers  by  Introducing  surface-dependent  nonlinear  current- 
potential  relationships  (polarization  curves),  as  shown  In  figure  9(c).  (It 
Is  anticipated  that  these  polarization  curves  may  vary  with  time.)  The 

current  flow  problem  across  the  electrode-electrolytic  Interface  appears  to 
be  conveniently  modeled  using  convection  elements  Instead  of  conduction 
elements.1*  Like  the  convection  boundary  problem  In  thermal  science,  the 
determining  factors  for  electrode  performance  are  su  face  area,  film 
resistance  location,  and  local  potential  difference.  Thi  boundary 

conditions  on  the  electrode  surfaces  are  constructed  by  assigning  the  anodic 
and  cathodic  potential  In  the  interfaces  to  the  uncoupled  open-circuit 

corrosion  potentials  of  the  constituents.  As  depicted  In  figure  9(c),  two 
nonlinear  film  elements  (CHBOY)  are  connected  to  the  electrolyte.  The 

bottom  film  elements  are  bounded  by  the  measured  open-circuit  potentials  for 
the  anodic  and  cathodic  surfaces  physically  described  In  figure  S.1**2^ 

An  additional  requirement  imposed  on  the  model  Is  that  the  divergence  of 
the  current  density  vector,  T,  must  be  equal  to  zero.  In  this 

investigation,  special  adjustments  were  made  to  the  COSNIC/NASTRAN  finite 
element  program  to  fulfill  this  condition.  The  COSMIC/NASI  RAN  program  and 

related  postprocessing  program  are  resident  on  a  VAX -11/780  computer  at  the 

Naval  Underwater  Systems  Center  In  New  London.  Connecticut. 

A  number  of  planes  were  Identified  In  order  to  more  easily  retrieve 
current  densities  and  voltage  Information  from  the  elements  and  grid  points, 
respectively.  Because  the  measurements  for  the  current  density  were 

performed  at  650  and  1150  i*i  above  the  electrode  surfaces.  It  was 
desirable  to  obtain  calculated  current  density  Information  at  the  centroids 
of  the  associated  elements.  The  lowest  plane  was  assumed  to  exist  at 

Z  *  0.  The  centroid  for  the  eleemnts  coincides  with  the  elevation  of  the 

measured  values  (Z  -  650  iae).  Similarly,  the  other  elements  from  the 

finite  element  model  have  a  centroldal  vertical  distance  of  Z  *  1125  yin, 

while  the  experimental  elevation  value  was  1150  ym.  The  topmost  plane 
(electrolyte/air  Interface)  Is  located  at  Z  »  9500  t*- 

It  was  determined  In  the  course  of  the  analysis  that  the  current  density 
of  the  6  direction  was  Indeed  zero  normal  to  the  R  Z  planes  for  o  =  0* 
and  e  =•  36*.  Thus,  the  boundary  conditions  were  satisfied.  Over  the 

entire  A6  slice,  the  current  density  was  identically  zero  because  the 
electrical  potentials  at  e  -  O’  and  b  -  36*  were  Identical.  However, 
within  the  finite  element  model,  very  small  current  densities  (approximately 
less  than  5  percent  of  Jj  and  Jp  in  the  O  direct ior)  were  evident. 
This  amount  of  error  is  expected  for  the  relatively  simple  interpolation 
functions  chosen.  Increased  numerical  resolution  can  be  achieved  by  finer 
discretization  in  areas  of  expected  high  gradient  variation  and  the  possible 
introduction  of  higher  order  Interpolation  functions  for  the  finite  elements. 
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CONCLUSIONS 


The  three-dimensional  current  density  structure  In  the  nearfield  of  an 
Iron-copper  galvanic  structure  was  characterized  by  finely  Measured  SVtl's 
at  Incremental  positions  In  an  electrolyte  In  space  at  Z  -  650  •««  and 
Z  *=  1125  above  the  electrode  surface.  A  corresponding  finite  element 
current  density  field  was  calculated  for  positions  In  space  similar  to  those 
that  were  measured.  Three-dliaenslonal  current  density  vectors 
J(JX,  Jy.  Jz)  ware  calculated,  measured,  and  subsequently  plotted 
(figures  11  through  15)  to  create  a  polnt-by-polnt  comparison  of  vector 
variables  In  electrolytic  space.  The  comparison,  component -versus - 
component,  was  In  excellent  agreement  for  this  relatively  simple  galvanic 
configuration. 

Based  on  the  completed  experimental  electric  field  scans  and  the 
corresponding  finite  element  field  predictions.  It  appears  that  the  finite 
element  technique  presents  a  strong  analytical  tool  for  calculating  the 
nearfleld  electric  Intensity  distributions  around  active  electrode  surfaces 
in  homogeneous  media.  The  finite  element  model  described  In  this  report 
made  reference  to  an  *a  priori*  measured  set  of  uncoupled  nonlinear 
polarization  curves  for  pure  Iron  and  pure  copper  (figure  8).  With  a  local 
tangent  slope  definition  of  the  potential  difference  at  a  given  position  of 
the  electrode  surface,  the  current  densities  were  calculated  everywhere 
within  the  electrolyte  and  on  the  boundaries. 

The  measurement  technique  described  has  Introduced  an  Innovative  method 
for  providing  a  careful  comparison  with  theory  In  the  very  nearfleld  of 
active  galvanic  cells.  Further  advances  could  lead  to  a  method  for 
predicting  the  onset  of  crevice  corrosion  and  pitting  In  field  situations. 
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Figure  1.  A  Schematic  of  the  Iron-Copper  Galvanic  Couple 
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Figure  6.  Schematic  of  the  Probe  Tip 
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Figure  7.  The  Geometric  Arrangement  of  the  Three-Dimensional  Scans 


Figure  8.  Potentiostatlc  Polarization  Curves  for  the  Iron  and 
Copper  Constituents  of  the  6alvanic  Couple 
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BOUNDARY  CONDITIONS  OF 
THE  THREE-DIMENSIONAL 
FINITE  ELEMENT  MESH  (36°) 

(a) 

NONLINEAR  SOLUTIONS  TO. 

WHERE 

dp  z  = 

V*(rV<I>  =  —  TT  O)  J  = 

-  0  T  <t>  = 

J  =  p  s 

t  = 

USING  MEASURED  CONSTITUENT  V  = 

POLARIZATION  CURVES 


CONDUCTIVITY  (S/m) 
CURRENT  DENSITY  (A/m2) 
POTENTIAL  (VOLTS) 
CHARGE  DENSITY  (C/m*) 
TIME  (s) 

DEL  OPERATOR 


(b) 


BOUNDARY  CONDITIONS  OF  MODELED 
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Figure  9.  Summary  of  the  Finite  Element  Model:  Geometric  Description, 
Fundamental  Electrical  Laws,  and  Boundary  Conditions 
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Figure  10.  Calculated  and  Measured  J2  Components  of  the  Ionic 
Current  Oensity  at  650-ym  Elevation 


Figure  11.  Calculated  and  Measured  Components  of  the  Ionic 
Current  Oensity  at  1125-ym  Elevation  (Measurement  at  1150  urn) 
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